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We explore perturbations about a Friedmann-Robertson- Walker background in Chern-Simons 
gravity. At large momenta one of the two circularly polarized tensor modes becomes ghostlike. We 
argue that nevertheless the theory does not exhibit classical runaway solutions, except possibly in 
the relativistic nonlinear regime. However, the ghost modes cause the vacuum state to be quantum 
mechanically unstable, with a decay rate that is naively infinite. The decay rate can be made finite 
only if one interprets the theory as an effective quantum field theory valid up to some momentum 
cutoff A, which violates Lorentz invariance. By demanding that the energy density in photons 
created by vacuum decay over the lifetime of the Universe not violate observational bounds, we 
derive strong constraints on the two dimensional parameter space of the theory, consisting of the 
cutoff A and the Chern-Simons mass. 



I. INTRODUCTION AND SUMMARY 

General relativity has held up well to various tests over 
the years from experiments and astronomical observa- 
tions [1], and is considered a pillar of standard cosmol- 
ogy. However, it is interesting to consider modifications 
to the theory, particularly in light of the observed accel- 
eration of the Universe [2j. One useful approach is to 
think of Einstein gravity as an effective field theory, and 
consider higher order corrections to the Einstein-Hilbcrt 
action, either involving the metric alone or involving an 
additional posited scalar field. The goal then becomes 
to calculate the corrections to general relativity arising 
from these higher order terms, and using experiments to 
set bounds on the couplings parameterizing them. 

One such extension to general relativity is Chern- 
Simons gravity [3J 0], where one assumes the existence 
of a canonical scalar field, coupled to gravity through a 
parity violating term. The theory is described by the 
action 



S — SeH + ScS + Stf + S'mat, 



(1.1) 



where the various terms are respectively the Einstein- 
Hilbert term 

(1.2a) 



the Chern-Simons term 
1 



Scs = / d 4 x^ 



-g -d *RR, 



(1.2b) 



the scalar term 
Sa = — — 



(fixyflj [V a tfV a t? + 2V{&)], (1.2c) 



and £ m at describes any other matter present. In these 
expressions m 2 = (87rG) _1 is the square of the reduced 
Planck mass, g is the determinant of the metric, R is the 
Ricci scalar and a a coupling constant with dimensions 



of inverse mass. Also V($) is an arbitrary potential and 
the Pontryagin density is defined as 



'RR 



1 



cdef no nb 

e it be ftt acd , 



(1.3) 



where e cde f is the four dimensional Levi-Civita tensor. 
For purposes of the present discussion, we will assume 
that matter is minimally coupled to the metric. 

Consider now the dynamics of Chern-Simons grav- 
ity in perturbation theory about a Friedman-Robertson- 
Walker (FRW) cosmological background. N ow in the 
limit $ = constant, the Chern-Simons term (1.2b) re- 



duces to a surface term and we recover Einstein's equa- 
tions for the space-time dynamics. Therefore, in the ef- 
fective field theory that describes the perturbations, the 
operators that arise from the Chern-Simons term must be 
suppressed by a mass scale that is related to the deriva- 
tive of the background scalar field. This mass scale is 
called the Chern-Simons mass m cs , and is defined by [5] 



(1.4) 



where the dot denotes a derivative with respect to time. 
General relativity is recovered in the limit m cs — > oo. Be- 
cause it is the Chern-Simons mass that enters into equa- 
tions describing observables, we choose to constrain it, 
rather than the more fundamental coupling a that ap- 
pears in the action (1.2b I. 



Past work constraining Chern-Simons gravity has uti- 
lized Solar System tests of general relativity. Measure- 
ment of Lense-Thirring precession by LAGEOS give the 
bound [5] 



> !CT 13 eV 



(1.5) 



A bound 10 11 times stronger has been claimed from bi- 
nary pulsar studies [5] , but the validity of this result has 
been questioned [7], and a corrected bound from binary 
pulsars is [7] 



> 5 x HT iU eV 



(1.6) 
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FIG. 1: Effective field theory parameter space consisting of 
the cutoff scale A and the Chern-Simons mass m cs , in eV. 
The lightly shaded region is excluded from binary pulsar ob- 
servations [7], and the darkly shaded region is excluded by 
consideration of vacuum decay. 



In this paper we study the vacuum stability of 
Friedmann- Robertson- Walker (FRW) solutions in Chern- 
Simons gravity as a function of the Chern-Simons mass 
parameter. In Section [TT] we consider tensor perturba- 
tions to the FRW metric. We show that for spatial mo- 
menta above the Chern-Simons mass scale, one of the 
two polarization modes is ghostlike and can decay to ra- 
diation. Requiring that the radiation produced over the 
lifetime of the Universe not exceed observational bounds 
will allow us to constraint the parameters of the theory. 

In Section |IV| we argue that the theory does not ex- 
hibit classical runaway solutions, despite the existence of 
ghostlike modes, except possibly in the relativistic non- 
linear regime. 



II. THE EXISTENCE OF GHOST GRAVITON 
MODES 

It is easily verified that an FRW metric, with scale 
factor and matter sources obeying the usual Friedmann 
equations, is a solution of the equations of motion of the 
theory (1.1) [5], as the FRW symmetries lead to a van- 
ishing Pontryagin density. Consider now linear perturba- 
tions about such a solution. As in general relativity, the 
symmetries of the background solution guarantee that 
perturbations can be decomposed into scalar, vector and 
tensor modes. 

In this paper we will focus on tensor perturbation 
modes, for which the perturbation to the scalar field van- 



ishes and the metric takes the form 



ds 2 = a 2 (77)[-d77 2 + (<% + hiAdx'dx?] 



(2.1) 



Here a(j]) is the scale factor, rj is conformal time and 
are comoving coordinates. We adopt the transverse and 
traceless gauge conditions, /i% = 0, dih^ — 0. Expand- 
ing the terms in the action ( 1.1 ) to quadratic order yields 
® 



drj d x 



0(h 3 ), 



(2.2) 



where e y is the Levi-Cevita symbol. We rewrite this 
action in terms of the Fourier transform of the metric 
perturbation, which is defined by 



hij(v,x) = / d 3 k hij(r},k) 



(2.3) 



where k is the comoving wavevector. We define the left 
and right circular polarization modes hAk(j]) by 



(2.4) 



A=L,R 



where n = k/k is a unit vector in the direction of prop- 
agation, k = \k\, and the polarization tensors e^-(n) 
satisfy the conditions 



e a ( e a ) 



(2.5b) 

-1 and Al = — 1. The action can now be 



25 



AB 



(2.5a) 



with Xji = 
written as 



q _ ^ 



f dr] d 3 k ° 2 (»7) 

X A k <x& >rj 
a 2 (r)) 772p 



[\h Ak , v \ 2 - k 2 \h <~ 



lAk 



(2.6) 



This can be recast using the Chern-Simons mass scale 
CI as 



S = T f [d v d 3 k £ a\n) 



A—L,R 



1 + A, 



%>hy 



[\h Ak , v \ 2 ~ k 2 \h Ak \ 2 ] , (2.7) 



wherefcphy = k/a is the physical wavenumber. The ac- 
tion (2.7) is the usual action for tensor perturbations in 



FRW, except for the momentum dependent correction 
factor (1 + A^fcphy/mcs). This factor becomes negative 



3 



for the left handed polarization modes when fc p hy > to cs , 
giving rise to a kinetic term with the wrong sign, i.e. a 
ghost mode. 



The action (2.7) can be simplified by changing the nor- 

;s to atta 

^ A fcphy s 



malization of the graviton modes to attain canonical nor- 
malization. We define 



s A (k) = sgn 1 



(2.8) 



which is +1 for normal modes and —1 for ghost modes. 
We define the mass scale 



1 



A fcphy s 



(2.9) 



and the canonically normalized graviton field modes 

K%% = mJi Ak . (2.10) 



The action (2.7) can now be written as 



5 = i f dr) d 3 k J2 a 2 (r!)e A {k) 



i, can 
n Ak,r) 



A=L,R 



can 
Ak 



.(2.11) 



Note that from a classical point of view the existence of 
ghost modes does not necessarily imply any inconsistency 
of the theory. The Hamiltonian of the theory may or 
may not be unbounded below; addressing this question 
would require a nonlinear analysis beyond the scope of 
this paper. Also it is not known at present whether the 
theory possesses a well posed initial value formulation: 
the sign flip at fc p hy S = to cs may be a hint that it does 
not. However, the ghost modes are a significant problem 
when quantum mechanical effects are taken into account, 
as we discuss in the next section. 

One might also expect that the ghost graviton modes 
would be associated with fluxes of negative energy. 
Rather surprisingly, this is not so; a computation of the 
effective stress energy tensor for gravitational waves at fu- 
ture null infinity shows that all the graviton modes carry 
positive energy [ID] . 



III. CONSTRAINTS FROM VACUUM DECAY 

We now specialize to perturbation modes today which 
are deep inside the horizon, that is, k 3> Hq, where Hq = 
a yrj /a 2 is the Hubble parameter. We also assume that 
to cs ^ > Hq. In this limit, we can neglect in the action 
(2.11) the time dependence of the prefactor a(rj) 2 , and 



also the term proportional to the time derivative of to*. 
The result is just the standard action for graviton modes 
in Minkowski spacctimc, except for the sign flip for the 
ghost modes 1 . 



1 We have omitted the action describing the interaction of the 
graviton field modes with matter fields. This interaction takes 



Because of this sign flip, decay of the vacuum of the 
theory is kinematically allowed. The vacuum decay rate 
per unit volume T is naively infinite, because of the in- 
finite phase space available for the decay products that 
arises from Lorentz invariance 2 . Rather than ruling out 
the theory outright because of the presence of g host 
modes, we adopt the viewpoint that the action (2.7) de- 



fines an effective quantum field theory with some effec- 
tive cutoff A on the physical wavenumber fc p h ys in cos- 
mological rest frame. Thus, our cutoff explicitly vio- 
lates Lorentz invariance; such a violation is inevitable 
if one wants to obtain a finite vacuum decay rate. Our 
viewpoint and treatment follow similar analyses of scalar 
ghost fields in cosmological models with equation of state 
parameter w that satisfies w < — 1 [TTJ|T2]. As in those 
analyses, we will find that stringent constraints on the 
parameters of the theory can be obtained by demanding 
that the total number of photons produced from vacuum 
decay over the lifetime of the Universe not be in conflict 
with observations. 

Our theory is now parameterized by two parameters, 
the cutoff A and the Chern-Simons mass to cs . Now since 
the ghost modes arise only for momenta fc p hy S satisfy- 
ing fcphys > to cs , there will be no ghost modes in the 
region m cs > A of parameter space. Therefore our argu- 
ments about vacuum decay do not constrain that region 
of parameter space; we focus from now on on the com- 
plementary region m cs < A (see Fig. [T]) . We note that 
many papers on Chern-Simons gravity implicitly work in 
this regime m cs < A, since they involve effects arising on 
scales fcphys ~ to cs , assumed to be within the domain of 
validity of the theory. 

If the vacuum perturbatively decays to stable parti- 
cles, we can estimate the number of particles produced 
over the lifetime of the Universe from such a process. 
Given observational bounds on the energy density in this 
particle species, we can then constrain the decay rate. 
The strongest bounds will therefore come from produc- 
tion of particles whose energy density has been reliably 
measured, and which can be produced at a low order in 
perturbation theory. We argue that the photon is the 
best such candidate. 

We next make an order of magnitude estimate of the 
decay rate T per unit volume to photons in the regime 
to cs <C A. The action for electromagnetism contains an 
interaction of the form 

S int ~ I d A xh(dA) 2 ~ — / d 4 xh can (OA) 2 , (3.1) 

J TO* J 

where A a is the 4-vector potential, h is the metric per- 



the standard form when written in terms of h^, but acquires 
correction factors of m p /m* when written in terms of h'TS- 
2 The full theory of the perturbations about FRW, including cou- 
pling to matter, does violate Lorentz invariance because of the 
fcphys dependence of the factor 1 + A a fc p h ys /m cs . However this 
violation does not affect the accessible volume of phase space. 
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FIG. 2: Spontaneous production of photons and ghost gravi- 
tons from vacuum 



turbation, and /i can is the canonically normalized version. 
Because the graviton is ghostlike, the process 



.977 



(3.2) 



is kinematically allowed, where g is a left polarized gravi- 
ton and 7 is a photon (Fig. [2| . Hence graviton ghosts can 
decay to photons at first order in perturbation theory. 

Next, the coefficient 1/m* in the interaction (3.1) de- 
pends on the wavenumber k. However, for the purposes 
of our order of magnitude estimate, it will be sufficient 
to evaluate this coefficient at fc p h ys ~ A, since most of 
the decays will be at /c p h ys ~ A. Therefore we can treat 
m* as a constant, 



A 



(3.3) 



where we have used m cs <C A and specialized to ghost 
modes. The decay rate Y per unit time per unit volume 
must be proportional to the square of the coefficient of 
the operator, so V oc 1/m 2 . The constant of proportion- 
ality in this relation must be some function of A, and 
from dimensional analysis it now follows that 



A 5 



(3.4) 



since T has dimension (mass) 4 in units with H = c = 1. 

Consider now the production of photons by this mech- 
anism over the lifetime of the Universe. We make the 
idealization that at any epoch, all the photons produced 
have an energy of exactly A. Then the evolution of the 
number n(fc, t) of photons per unit logarithmic wavenum- 
ber k per proper volume is given by 



= TA8(k/a — A), 



(3.5) 



where a is the scale factor with a = 1 today. Solving this 
equation gives for the spectrum today 



n(k) = 9(A - k) 



(3.6) 



where is the step function, H = a/ a is the Hubble 
parameter, and <z*(fe) and are the values of a and 

H at k/a = A. The result (3.6) is easy to understand: the 



factor of is the length of time during which photons 
of present-day energy ~ k are produced, and the factor 
of al is the volume expansion factor since then. For a 
ACDM cosmology with H 2 = H%(Q M /a 3 + 1 - fl M ) the 
spectrum becomes 



i(fc) = e(A-fc) 



M 



M 



A 



1-1/2 



(3.7) 



It follows that the energy density per unit logarithmic 
wavenumber dE/d 3 xd\nk ~ kn(k) is peaked at k ~ A 
with a peak value of 



dE TA 
d 3 xd In k H 



A 6 



(3.8) 



We now compare the prediction (3.8) with observa- 
tional data. Observations in various wavelength bands 
from radio to gamma rays gives the upper bound 3 on the 
background radiation spectrum 



dE 



d?xd\nk 



< 



(2.5 x l(T 5 eV) 4 



(3.9) 



for lCT 10 eV < k < lOGeV [HCS]. Combining this with 
the prediction (3.8) yields the constraint 



m cs A 6 < (3.3 eV) 7 , (3.10) 

for 10 -10 eV < A < lOGeV, where we have used m p — 
2.4 x 10 27 eV and H = 1.5 x 1(T 33 eV. Figure [l] shows 



the constraints (1.6) and (3.10) on the (m cs , A) parameter 



space. It can be seen that the combined constraints rule 
out a large part of the region in parameter space with A > 
m cs in which interesting Chern-Simons phenomenology 
occurs, except for the window 10 -10 eV < m cs < 3eV. 

One can also compute the constraints obtained from 
vacuum decay to gravitons, in which one ghost graviton 
and two normal gravitons are produced. The interaction 
Lagrangian is 



£mt ~ m 2 h(dh) 



p -h (dh ) 2 



(3.11) 



The analysis now proceeds as before, with the resulting 
energy density in gravitons given by 



1 



dE 



A 4 ? 



gw 



iti 2 Hq d 3 xd\nk m^H^ 



(3.12) 



at k ~ A. Now observations of the cosmic mi- 
crowave background give an integrated upper bound 
/ dlnfcO gw < 10~ 5 for k > lQ- 30 eV [T5]. However this 



3 This upper bound is exceeded by about two orders of magnitude 
by the CMB, but the CMB is known to be thermal to about 1 
part in 10 5 , so the bound effectively holds. 
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bound applies only to the gravitons produced before re- 
combination, whose energy density will be smaller than 



the estimate (3.12) by a factor of the ratio of the age 



of the Universe at recombination to the age of the Uni- 
verse today. This factor is approximately z rc ^ 2 , where 
z lec ~ 1000 is the redshift of recombination. We therefore 
obtain the constraint 

m 3 s A 4 < z^m^H 3 (10" 5 ) ~ (32 eV) 7 (3.13) 

for A > 10~ 30 eV, which is slightly weaker than the con- 
straint (3.10) from photons. 



IV. CLASSICAL RUNAWAY SOLUTIONS 

We next turn to a different question, the possible ex- 
istence of classical runaway solutions that can occur in 
a theory whose Hamiltonian is unbounded below. We 
will argue that such solutions do not arise in Chern- 
Simons gravity except possibly in the nonlinear relativis- 
tic regime h ~ 1 where our analysis is invalid anyway. 

Generically, a Hamiltonian which is unbounded from 
below will only exhibit runaway solutions in the regime 
where the interaction terms in the Lagrangian are com- 
parable to the negative kinetic energy terms [TT] . We use 
this criterion to estimate the required occupation number 
of a ghost graviton mode for such a runaway behavior to 
develop in our system. 

No runaway solutions will plague the effective field 



theory (2.7) because the ghost and non-ghost fields are 
decoupled to linear order. However, interaction terms 
will appear if we expand the action flXTXj ) to higher than 
quadratic order. Expanding up to quartic order we ex- 
pect to find interaction terms of the form 



So 



k 2 hjh, 



S4 ~ J d i x m 2 ^k 2 h 2 L h 2 R 



(4.1a) 



(4.1b) 



Consider a single localized wave packet mode with char- 
acteristic size A. We would like to estimate the number 
of quanta N for an instability to form. The energy from 
the kinetic term is 



E K = [ d 3 x ml h 2 ~ \m 2 h 2 ~ — . 
2 J p p A 

Therefore the field fluctuation is approximately 
. N 1 ' 2 



(4.2) 



(4.3) 



Each of the interaction terms will therefore contribute 
terms to the energy on the order of 

r jv 3 / 2 

£3 ~ / d 3 x m 2 k 2 h 3 ~ , (4.4a) 



N 2 

d 3 xm 2 k 2 h 4 ~^-^. (4.4b) 
1 \ 6 m± 



Requiring that the kinetic energy be less than the inter- 
action energy 



E3 > E K , £4 > E K , 

and taking the wavelength A ~ m^ 1 yields 



N > 



» 1. 



We deduce from (4.3) that in this regime 

h > 1. 



(4.5) 



(4.6) 



(4.7) 



We conclude the ghost modes generated by perturbations 
to FRW do not exhibit runaway solutions within the do- 
main of validity of our analysis. 



V. OUTLOOK 

We have shown that if the effective field theory cutoff 
A is above the Chern-Simons mass scale m cs , the vac- 
uum of Chern-Simons gravity is unstable and can de- 
cay to photons. Our estimate of the photon production 
rate was used to set constraints on the parameter space. 
One might imagine improving this bound by being less 
cavalier in the estimation of the decay rate and num- 
ber density. However, because the production rate goes 
like a high power of the cutoff scale, orders of magni- 
tude difference in the number density yield very modest 
improvements in final constraints on A and m cs . 

Whether the full nonlinear theory has a Hamiltonian 
which is unbounded below is an interesting open ques- 
tion. 

Finally, we note that it is possible that the decay of 
the vacuum could be a lot faster than indicated by our 
perturbative calculations, if it is mediated by a non- 
perturbative process. One might imagine a tunneling 
process similar to the decay of false vacua in scalar field 
theories via bubble nucleation. In particular, a homo- 
geneous unstable region with small Chern-Simons might 
be expected to produce a stable endstate consisting of 
a homogeneous region with a Chern-Simons mass larger 
than the cutoff. However any such transition could not 
proceed via a spherically symmetric process, since the 
Chern-Simons term does not contribute to the dynamics 
in spherical symmetry. 
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